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Here are treated the elements of the theory of the motion of fluids in general, the whole matter being reduced to 
this: given a mass of fluid, either free or confined in vessels, upon which an arbitrary motion is impressed, and 
which in turn is acted upon by arbitrary forces, to determine the motion carrying forward each particle, and at 
the same time to ascertain the pressure exerted by each part, acting on it as well as on the sides of the vessel. 
At first in this memoir, before undertaking the investigation of these effects of the forces, the Most Famous 
Authofl carefully evaluates all the possible motions which can actually take place in the fluid. Indeed, even 
if the individual particles of the fluid are free from each other, motions in which the particles interpenetrate 
are nevertheless excluded, since we are dealing with fluids that do not permit any compression into a narrower 
volume. Thus it is clear that an arbitrary small portion of fluid cannot receive a motion other than the one which 
constantly conserves the same volume; even though meanwhile the shape is changed in any way. It would hold 
indeed, as long as no elementary portion would be compressed at any time into a smaller volume; furthermore 
if the portion expanded into a larger volume, the continuity of the particles was violated, these were dispersed 
and no longer clinged together, such a motion would no longer pertain to the science of the motion of fluids; 
but individual droplets would separately perform their motion. Therefore, this case being excluded, the motion 
of the fluids must be restricted by this rule that each small portion must retain for ever the same volume; and 
this principle restricts the general expressions of motion for elements of the fluid. Plainly, considering an 
arbitrary small portion of the fluid, its individual points have to be carried by such a motion that, when at a 
moment of time they arrive at the next location, till then they occupy a volume equal to the previous one; thus 
if, as usual, the motion of a point is decomposed parallel to fixed orthogonal directions, it is necessary that 
a certain established relation hold between these three velocities, which the Author has determined in the first part. 

In the second part the author proceeds to the determination of the motion of a fluid produced by arbitrary forces, 
in which matter the whole investigation reduces to this that the pressure with which the parts of the fluid at each 
point act upon one another shall be ascertained; which pressure is denoted most conveniently, as customary for 
water, by a certain height; this is to be understood thus, that the each element of the fluid sustains a pressure the 
same as if were pressed by a heavy column of the same fluid, whose height is equal to that amount. Thus, in such 
way in each point of the fluid the height referring to the state of the pressure will be given; since it is not equal 
to the one in the neighbourhood, it will perturb the motion of the elements. But this pressure depends as well on 
the forces acting on each element of the fluid, as on those, acting in the whole mass; thus, by the given forces, the 
pressure in each point and thereupon the acceleration of each element - or its retardation - can be assigned for the 
motion, all which determinations are expressed by the author through differential formulae. But, in fact, the full 
development of these formulas mostly involves the greatest difficulties. But nevertheless this whole theory has 
been reduced to pure analysis, and what remains to be completed in it depends solely upon subsequent progress 
in Analysis. Thus it is far from true that purely analytic researches are of no use in applied mathematics; rather, 
important additions in pure analysis are now required. 



I. FIRST PART 

1. Since liquid substances differ from solid ones by the 
fact that their particles are mutually independent of each other, 



"This is an English adaptation by Walter Pauls of Euler's 
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they can also receive most diverse motions; the motion per- 
formed by an arbitrary particle of the fluid is not determined 
by the motion of the remaining particles to the point that it 
cannot move in any other way. The matter is very different 
in solid bodies, which, if they were inflexible, would not un- 
dergo any change in their shape; in whatsoever way they be 
moved, each of their particles would constantly keep the same 
location and distance with respect to other particles; it thus 
follows that, the motion of two or, if necessary, three of all the 
particles being known, the motion of any other particle can be 
defined; furthermore the motion of two or three particles of 
such a body cannot be chosen at will, but must be constrained 
in such a way that these particles preserve constantly their po- 
sitions with respect to each other. 1 

2. But if, moreover, solid bodies are flexible, the motion 
of each particle is less constrained: because of the bending, 



Here Euler refers to the motion of rigid solid bodies treated previously in 
Euler, 1750. 
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the distance as well as the relative position of each particle 
can be subject to changes. However, the manner itself of the 
bending constitutes a certain law which various particles of 
such a body have to obey in their motion: certainly what has 
to be taken care of is that the parts that experience in their 
neighborhood such a strong bending with respect to each other 
are neither torn apart from the inside nor penetrate into each 
other. Indeed, as we shall see, impenetrability is demanded 
for all bodies. 

3. In fluid bodies, whose particles are united among them- 
selves by no bond, the motion of each particle is much less 
restricted: the motion of the remaining particles is not de- 
termined from the motion of any number of particles. Even 
knowing the motion of one hundred particles, the future mo- 
tion permitted to the remaining particles still can vary in in- 
finitely many ways. From which it is seen that the motion of 
these fluid particles plainly does not depend on the motion of 
the remaining ones, unless it be enclosed by these so that it is 
constrained to follow them. 

4. However, it cannot happen that the motion of all parti- 
cles of the fluid suffers no restrictions at all. Furthermore, one 
cannot at will invent a motion that is conceived to occur for 
each particle. Since, indeed, the particles are impenetrable, 
it is immediately clear that a motion cannot be maintained in 
which some particles go through other particles and, accord- 
ingly, penetrate each other: also, because of this reason such 
motion certainly cannot be conceived to occur in the fluid. 
Therefore, infinitely many motions must be excluded; after 
their determination the remaining ones are grouped together. 
It is seen worthwhile to define them more accurately regard- 
ing the property which distinguishes them from the previous 
ones. 

5. But before the motion by which the fluid is agitated at 
any place can be defined, it is necessary to see how every mo- 
tion, which can definitely be maintained in this fluid, be rec- 
ognized: these motions, here, I will call possible, which I will 
distinguish from impossible motions which certainly cannot 
take place. We must then find what characteristic is appro- 
priate to possible motions, separating them from impossible 
ones. When this is done, we shall have to determine which 
one of all possible motions in a certain case ought actually 
to occur. Plainly we must then turn to the forces which act 
upon the water, so that the motion appropriate to them may be 
determined from the principles of mechanics. 

6. Thus, I decided to inquire into the character of the possi- 
ble motions, such that no violation of the impenetrability can 
occur in the fluid. I shall assume the fluid to be such as never 
to permit itself to be forced into a lesser space, nor should its 
continuity be interrupted. Once the theory of fluids has been 
adjusted to fluids of this nature, it will not be difficult to ex- 
tend it also to those fluids whose density is variable and which 
do not necessarily require continuity. 2 

7. If, thus, we consider an arbitrary portion in such a fluid, 



the motion, by which each of its particles is carried has to be 
set up so that at each time they occupy an equal volume. When 
this occurs in separate portions, any expansion into a larger 
volume, or compression into a smaller volume is prohibited. 
And, if we turn attention to this only property, we can have 
only such motion that the fluid is not permitted to expand or 
compress. Furthermore, what is said here about arbitrary por- 
tions of the fluid, has to be understood for each of its elements; 
so that the volume of its elements must constantly preserve its 
value. 

8. Thus, assuming that this condition holds, let an arbi- 
trary motion be considered to occur at each point of the fluid; 
moreover, given any element of the fluid, consider the brief 
translations of each of its boundaries. In this manner the vol- 
ume, in which the element is contained after a very short time, 
becomes known. From there on, this volume is posed to be 
equal to the one occupied previously, and this equation will 
prescribe the calculation of the motion, in so far as it will be 
possible. Since all elements occupy the same volumes during 
all periods of time, no compression of the fluid, nor expansion 
can occur; and the motion is arranged in such a way that this 
becomes possible. 

9. Since we consider not only the velocity 3 of the motion 
occuring at every points of the fluid but also its direction, both 
aspects are most conveniently handled, if the motion of each 
point is decomposed along fixed directions. Moreover, this de- 
composition is usually carried out with respect to two or three 
directions: 4 the former is appropriate for the decomposition, 
if the motion of all points is completed in the same plane; but 
if their motion is not contained in the same plane, it is ap- 
propriate to decompose the motion following three fixed axes. 
Because the latter case is more difficult to treat, it is more con- 
venient to begin the investigation of possible motions with the 
former case; once this has been done, the latter case will be 
easily completed. 

10. First I will assign to the fluid two dimensions in such 
a way that all of its particles are now not only found with 
certainty in the same plane, but also their motion is performed 
in it. Let this plane be represented in the plane of the table 
(Fig. 1), let an arbitrary point I of the fluid be considered, 
its position being denoted by orthogonal coordinates AL = 
x and LI = y. The motion is decomposed following these 
directions, giving a velocity Im — u parallel to the axis AL 
and In — v parallel to the other axis AB: so that the true future 
velocity of this point is = \J\uu + vv), and its direction with 
respect to the axis AL is inclined by an angle with the tangent 

v_ 
u ' 

11. Since the state of the motion, presented in a way which 
suits the each point of the fluid, is supposed to evolve, the 
velocities u and v will depend on the position I of the point 
and will therefore be functions of the coordinates x and y. 



2 See the English translation of "General laws of the motion of fluids" in 
these Proceedings. 



3 Meaning here the absolute value of the velocity. 

4 Depending on the dimension: Euler treats both the two- and the three- 
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Thus, we put upon a differentiation 

du = Ldx + Idy and dv = Mdx + mdy, 

which differential formulas, since they are complete, 5 satisfy 
furthermore ^ = ^ and 7^ = 757- Here it is to note that 
in such expression 4J^, the differential of L itself or dh, is 
understood to be obtained from the variability with respect 
to y; in similar manner in the expression dl/dx, for dl the 
differential of I itself has to be taken, which arises if we take 
x to vary. 

12. Thus, it is in order to be cautious and not to take in 
such fractional expressions 4^, 4^, 4r^, and 4^ the numera- 

r ay 7 ax 7 ay 7 ax 

tors dL, dl, dM, and dm as denoting the complete differentials 
of the functions L, I, M and m; but constantly they designate 
such differentials that are obtained from the variation of only 
one coordinate, obviously the one, whose differential is repre- 
sented in the denominator; thus, such expressions will always 
represent finite and well defined quantities. Furthermore, in 
the same way are understood L= 4^, I = 4^, M= ^- and 

J ax 7 ay 7 ax 

m = 4^; which notation of ratios has been used for the first 

dy 7 

time by the most enlightened Fontaine, 6 and I will also apply 
it here, since it gives a non negligible advantage of calculation. 

13. Since du = Ldx+ldy and dv = Mdx+mdy, here it is 
appropriate to assign a pair of velocities to the point which is 
at an infinitely small distance from the point I; if the distance 
of such a point from the point I parallel to the axis AL is dx, 
and parallel to the axis AB is dy, then the velocity of this point 



5 Exact differentials. 

6 A paper "Sur le calcul integral" containing the notation & for the partial 
derivative of / with respect to x was presented by Alexis Fontaine des 
Bertins to the Paris Academy of Sciences in 1738, but it was published 
only a quarter of a century later (Fontaine, 1764). Nevertheless, Fontaine's 
paper was widely known among mathematicians from the beginning of 
the 1740s, and, particularly, was discussed in the correspondence between 
Euler, Daniel Bernoulli and Clairaut; cf. Euler, 1980: 65-246. 



parallel to the axis AL will be u + Ldx + Idy; furthermore, the 
velocity parallel to the other axis AB is v + Mdx + mdy. Thus, 
during the infinitely short time dt this point will be carried 
parallel to the direction of the axis AL the distance dt(u + 
Ldx + Idy) and parallel to the direction of the other axis AB 
the distance dt(v + Mdx + mdy). 

14. Having noted these things, let us consider a triangular 
element Imn of water, and let us seek the location into which 
it is carried by the motion during the time dt. Let Im be the 
side parallel to the axis AL and let In be the side parallel to 
the axis AB: let us also put Im = dx and In = dy; or let the 
coordinates of the point mbe x + dx and y; the coordinates of 
the point n be x and y + dy. It is plain, since we do not define 
the relation between the differentials dx and dy, which can be 
taken negative as well as positive, that in thought the whole 
mass of fluid may be divided into elements of this sort, so that 
what we determine for one in general will extend equally to 
all. 

15. To find out how far the element Imn is carried during 
the time dt due to the local motion, we search for the points 
p, q and r, to which its vertices, or the points I, m and n are 
transferred during the time dt. Since 





of point I 


of point m 


of point n 


Velocity w.r.t. AL= 


u 


u + Ldx 


u + Idy 


Velocity w.r.t. AB= 


V 


v + Mdx 


v + mdy 



in the time dt the point I reaches the point p, chosen such that: 

AP - AL = udt and Pp - LI = vdt. 
Furthermore, the point m reaches the point q, such that 
AQ-AM = (u + Ldx)dt and Qg-Mm = (v + Mdx)dt. 
Also, the point n is carried to r, chosen such that 

AR — AL = (u + ldy)dt and Rr — Ln = (v + mdy)dt. 

16. Since the points /, m and n are carried to the points p, q 
and r, the triangle Imn made of the joined straight lines pq, pr 
and qr, is thought to be arriving at the location defined by the 
triangle pqr. Because the triangle Imn is infinitely small, its 
sides cannot receive any curvature from the motion, and there- 
fore, after having performed the translation of the element of 
water Imn in the time dt, it will conserve the straight and 
triangular form. Since this element Imn must not be either 
extended to a larger volume, nor compressed into a smaller 
one, the motion should be arranged so that the volume of the 
triangle pqr is rendered to be equal to the area of the triangle 
Imn. 

17. The area of the triangle Imn, being rectangular at I, is 
= \dxdy, value to which the area of the triangle pqr should 
be put equal. To find this area, the pair of coordinates of the 
points p, q and r must be examined, which are: 

AP = x + udt; AQ = x + dx + (u + Ldx)dt; 

AR = x + (u + ldy)dt; Pp = y + vdt 

Qq = y + (v + Mdx)dt, Rr = y + dy + (v + mdy)dt 
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Then, indeed, the area of the triangle pqr is found from the 
area of the succeeding trapezoids, so that 

pqr — PprR + RrgQ — PpqQ. 

Since these trapezoids have a pair of sides parallel and per- 
pendicular to the base AQ, their areas are easily found. 

18. Plainly, these areas are given by the expressions 

PprR = ipR(Pp + Rr) 
RrqQ = ^RQ(Rr + Qq) 

Pp<zQ = ^PQ(Pp + Qg) 

By putting these together we find: 

Apqr = ipQ.Rr - ^RQ.Pp - ^PR.Qg 
Let us set for brevity 

AQ = AP + Q; AR = AP + R; Qg = Pp + q; and 
Rr = Pp + r, 

so that PQ = Q, PR = R, and RQ = Q - R, and we 
have Apqr = ±Q(Pp + r) - |(Q - R)Pp - ±R{Pp + q) or 
Apqr = iQ.r - \R.q. 

19. Truly, from the values of the coordinates represented 
before it follows that 

Q = dx + Ldxdt; q = Mdxdt 
R = Idydt; r = dy + mdydt, 

upon the substitution of these values, the area of the triangle 
is obtained 

pqr = -dxdy(l + Ldt)(l + mdt) — -Ml dxdydt 2 , or 

pqr = ^dxdy(l + Ldt + mdt + Lmdt 2 - Mldt 2 ). 

This should be equal to the area of the triangle Iran, that is 
= \dxdy; hence we obtain the following equation 

Ldt + mdt + Lmdt 2 - Mldt 2 = or 
L + to + Lmdt - Mldt = 0. 

20. Since the terms Lmdt and Mldt vanish for finite L and 
to, we will have the equation L + m = 0. Hence, for the mo- 
tion to be possible, the velocities u and v of any point I have 
to be arranged such that after calculating their differentials 

du = Ldx + Idy, and dv = Mdx + mdy, 

one has L + m = 0. Or, since L= ^ and m = 4^, the 

' ax dy 7 

velocities u and v, which are considered to occur at the point 
I parallel to the axes AL and AB, must be functions of the 
coordinates x and y such that ^ + ^ = 0, and thus, the 

criterion of possible motions consists in this that + 3^ = 



0; 7 and unless this condition holds, the motion of the fluid 
cannot take place. 

21. We shall proceed identically when the motion of the 
fluid is not confined to the same plane. Let us assume, to in- 
vestigate this question in the broadest sense, that all particles 
of the fluid are agitated among themselves by an arbitrary mo- 
tion, with the only law to be respected that neither condensa- 
tion nor expansion of the parts occurs anywhere: in the same 
way, we seek which condition should apply to the velocities 
that are considered to occur at every point, so that the motion 
is possible: or, which amounts to the same, all motions that 
are opposed to these conditions should be eliminated from the 
possible ones, this being the criterion of possible motions. 

22. Let us consider an arbitrary point of the fluid A. To 
represent its location we use three fixed axes AL, AB and AC 
orthogonal to each other (Fig. 2). Let the triple coordinates 
parallel to these axes be AL= x, LI = y and = z; which 
are obtained if firstly a perpendicular \l is dropped from the 
point A to the plane determined by the two axes AL and AB; 
and then a perpendicular IL is drawn from the point I to the 
axis AL. In this manner the location of the point A is expressed 
through three such coordinates in the most general way and 
can be adapted to all points of the fluid. 




23. Whatever the later motion of the point A, it can be 
resolved following the three directions A/i, \v, Ao, parallel to 
the axes AL, AB and AC. For the motion of the point A we 
set 

the velocity parallel to the direction A/i = u, 
the velocity parallel to the direction \v = v, 
the velocity parallel to the direction Ao = w. 

Since these velocities can vary in an arbitrary manner for dif- 
ferent locations of the point A, they will have to be considered 
as functions of the three coordinates x, y and z. After differ- 



7 This is the two-dimensional incompressibility condition, which in a slightly 
different form has already been established by d'Alembert, 1752; cf. also 
Darrigol and Frisch, 2008:§ III. 
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entiating them, let us put to proceed 

du = Ldx + Idy + Xdz 
dv = Mdx + mdy + \xdz 
dw = Ndx + ndy + vdz. 

Henceforth the quantities L, I, A, M, m, p, N, n, v will be 
functions of the coordinates x, y and z. 

24. Because these formulas are complete differentials, we 
obtain as above 



dL 


dl 


dL 


dX 


dl 


d\ 


dy 


dx ' 


dz 


dx ' 


dz 


dy 


dM 


dm 


dM 


dp 


dm 


dp 


dy 


dx ' 


dz 


dx ' 


dz 


dy 


dN 


dn 


dN 


dv 


dn 


dv 


dy 


dx ' 


dz 


dx ' 


dz 


dy 



where it is assumed that the only varying coordinate is that 
whose differential appears in the denominator varies. 8 

25. Thus, this point A will be moved in the time dt by 
this threefold motion, which is considered to take place at the 
point X; hence it moves 

parallel to the axis AL the distance = udt 
parallel to the axis AB the distance = vdt 
parallel to the axis AC the distance = wdt 

The true velocity of the point A, denoted by = V, which 
clearly arises from the composition of this triple motion, is 
given in view of orthogonality of the three directions by 
V = \f(uu + vv + ww) and the elementary distance, which 
is traveled in the time dt through its motion, will be Vdt. 

26. Let us consider an arbitrary solid element of the fluid to 
see whereto it is carried during the time dt; since it amounts 
to the same, let us assign a quite arbitrary shape to that ele- 
ment, but of a kind such that the entire mass of the fluid can 
be divided into such elements; to investigate by calculation, 
let the shape be a right triangular pyramid, bounded by four 
solid angles A, p, v and o, so that for each one there are three 
coordinates 





of point A 


of point p 


of point v 


of point o 


w.r.t. AL 


X 


x + dx 


X 


X 


w.r.t. AB 


y 


y 


y + dy 


y 


w.r.t. AC 


z 


z 


z 


z + dz 



Since the base of this pyramid is Xpv = Imn = ^dxdy and 
the hight Ao = dz, its volume will be = ^dxdydz. 

27. Let us investigate, whereto these vertices A, p, v and 
o are carried during the time dt: for which purpose their three 
velocities parallel to the directions of the three axes must be 



The partial differential notation was so new that Euler had to remind the 
reader of its definition. 



considered. The differential values of the velocities u, v and 
w are given by 

Velocity of point A of point p of point v of point o 
w.r.t. AL u u + Ldx u + Idy u + Xdz 
w.r.t. AB v v + Mdx v + mdy v + pdz 
w.r.t. AC w w + Ndx w + ndy w + odz 

28. If we let the points A, p, v and o be transferred to the 
points 7r, $, p and a in the time dt, and establish the three 
coordinates of these points parallel to the axes, the small dis- 
placement parallel to these axes will be 

AP - A L = udt 

AQ-A M = (u + Ldx) dt 

AR-AL = (u + l dy)dt 

AS -A L = (u + A dz) dt 
Pp — LI = vdt 
Qq - Mm = (v + Mdx) dt 
Rr — L n = (v + mdy) dt 
Ss — LI = (v + p, dz) dt 
pir — I X = wdt 
q<& - mpi, = (w + N dx) dt 
r p — nv = (w + n dy) dt 
s a — I o = (w + v dz) dt 

Thus the three coordinates for these four points n, $, p and a 
will be 



AP = x + udt; 

PTT = Z + Wdt 

RQ = x + dx + (u + Ldx)dt; 
q$ = z+(w + Ndx)dt 
AR = i + (u + ldy)dt; 
rp — z + (w + ndy)dt 
AS = x + (u + Xdz)dt; 
sa = z + dz + (w + vdz)dt 



Pp = y + vdt; 

Qq = y+(v + Mdx)dt; 

Rr = y + dy + (v + mdy)dt; 

Ss = y + (v + pdz)dt; 



29. Since after the time dt has elapsed the vertices A, p, 
v and o of the pyramid are transferred to the points ir, $, p 
and a, ■n^pa defines a similar triangular pyramid. Due to the 
nature of the fluid the volume of the pyramid irQpa should 
be equal to the volume of the pyramid Xpvo put forward, that 
is = ^dxdydz. Thus, the whole matter is reduced to deter- 
mining the volume of the pyramid ir^pa. Clearly, it remains 
a pyramid, if the solid pqrir^pa is removed from the solid 
pqrir$pa; the latter solid is a prism orthogonally incident to 
the triangular basis pqr, and cut by the upper oblique section 

30. The other solid pqrirQpa can be divided by similarly 
into three prisms truncated in this manner, namely 

I. pqrsn^cr; II. prsTrpa; III. qrs$pa 

This has to be accomplished in such a way that 

—dxdydz — pqsn^a + prsirpcr + qrs^pu — pqrirQp. 
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Since such a prism is orthogonally incident to its lower base, 
and furthermore has three unequal heights, its volume is found 
by multiplying the base by one third of the sum of these 
heights. 

31. Thus, the volumes of these truncated prisms will be 

pq s 7T $(T = \-p q s (p it + q $ + s a ) 
o 

1 

pr s npa — -pr s (pir + r p + s a) 

O 

qr s $p a = -q r s (q $ + r p + s a) 

pqr 7r$p = \-p q r (pir + q $ + r p). 

Since pqr — pqs+prs + qrs, the sum of the first three prisms 
will definitely be small, or 

1 , , , 1 1 ^ 1 1 
-dxdydz = pir. qrs q<P.prs rp.pqs H — sa.pqr, 

6 3 3 3 3 

or 

dxdydz = Ipqr.sa — 2pqs.rp — 2prs.q& — 2qrs.pir. 

32. Thus, it remains to define the bases of these prisms: 
but before we do this, let us put 





AP + Q; 


Qq 


= Pp 


+ q; 


<?$ 


= pir 




AR = 


AP + R; 


Rr 


= Pp 


+ r ; 


rp 


= pir 


+ P 


AS = 


AP + S; 


Ss 


= Pp 


+ s; 


sa 


= PTT 


+ a 



in order to shorten the following calculations. After the substi- 
tution of these values, the terms containing pir will annihilate 
each other, and we shall have 

dxdydz = 2pqr.a — 2pqs.p — 2prs.<& 

so that the value of the bases to be investigated is smaller. 

33. Furthermore the triangle pqr is obtained by removing 
the trapezoid PpqQ from the figure PprqQ, the latter being the 
sum of the trapezoids PprR and RrqQ; from which it follows 
that 

Apqr = ipR(Pp+Rr) + iRQ(Rr+Q g )-ipQ(Pp+Qg); 

or, because of PR = R; RQ = Q — R; and PQ = Qwe shall 
have 

Apqr = ^R(Pp - Qq) + ^Q(Rr - Pp) = ^Qr - ^Rg. 
In the same manner we have 



And finally 
1 



-PR(Pp + Rr) + -RS(Rr 



Aprs 
or 

Aprs 

from where it follows that 



Ss) - -PS(Pp + Ss) 



^R(Pp + Rr) + -(S - R)(Rr + Ss) - ^S(Pp + Ss) 



Aprs = ^R(Pp - Ss) + is(Rr - Pp) = ^Sr - ^Rs. 

34. After the substitution of these values we will obtain 

dxdydz = (Qr - Rq)a + (Sq - Qs)p + (Rs - Sr)$; 
thus the volume of the pyramid ir&pa will be 

i(Qr - Rq)a + ^(Sq - Qs)p + i(Rs - Sr)$. 

From the values of the coordinates presented above in §. 28 
follows 

Q = dx + hdxdt q = Mdxdt $ = Ndxdt 
R = Idydt r = dy + m,dydt p = ndydt 
S = Xdzdt s = pdzdt a = dz + vdzdt. 

35. Since here we have 

Qr - Rq = dxdy{\ + Ldt + mdt + Lmdt 2 - Mldt 2 ) 
Sq -Qs = dxdz(-pdt-Lpdt 2 + MXdt 2 ) 
Rs - Sr = dydz(-\dt - mXdt 2 + Ipdt 2 ) 

the volume of the pyramid ir$pa is found to be expressed as 

' 1 +L dt +Lmdt 2 +Pmv dt 3 * 

+mdt -Ml dt 2 -Mlu dt 3 

+v dt +Lu dt 2 —Lnp dt 3 

+mudt 2 +MnXdt 3 

-np dt 2 -NmXdt 3 

-NA dt 2 +Nlp dt 3 J 



i dxdydz < 



> > 



which (volume), since it must be equal to that of the pyramid 
Xpvo — ^dxdydz, will satisfy, after performing a division by 
dt the following equation 9 

=L + m + v + dt(Lm + Lu + mu- Ml - NA - np) 
+ dt 2 (Pmv + MnX + Nip, - Lnp - Mlu - Nip). 



Apqs = -PS(Pp + Ss) H — SQ(Ss + Qg) PQ(Pp + Qq), 36. Discarding the infinitely small terms, we get this equa- 

"' tion: 10 L+m + u = 0, through which is determined the re- 

or lation between u, v and w, so that the motion of the fluid be 

Apqs = is(Pp + Ss) + 1(Q - S)(Ss + Qq) - ^Q(Pp + Qq), 



from where it follows that: 
1 



1. 



= -S(Pp - Qq) + ^Q(Ss - Pp) = ^Qs - -Sq. 



9 This is the calculation to which Euler refers in his later French memoir 
Euler, 1755. 

10 This is the three-dimensional incompressibility condition. 
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possible. Since L= m — ^ and = ^j, at an arbitrary 
point of the fluid A, whose position is defined by the three 
coordinates x, y and z, and the velocities u, v and w are as- 
signed in the same manner to be directed along these same 
coordinates, the criterion of possible motions is such that 

du dv dw 
dx dy dz 

This condition expresses that through the motion no part of the 
fluid is carried into a greater or or lesser space, but perpetually 
the continuity of the fluid as well as the identical density is 
conserved. 

37. This property is to be interpreted further so that at the 
same instant it is extended to all points of the fluid: of course, 
the three velocities of all the points must be such functions of 
the three coordinates x, y and z that we have ^| + ^ + = 
0: in this way the nature of those functions defines the motion 
of every point of the fluid at a given instant. At another time 
the motion of the same points may be howsoever different, 
provided that at an arbitrary point of time the property holds 
for the whole fluid. Up to now I have considered the time 
simply as a constant quantity. 

38. If however, we also wish to consider the time as vari- 
able so that the motion of the point A whose position is given 
by the three coordinates AL= x,Ll — y and IX = z has to 
be defined after the elapsed time t, it is certain that the three 
velocities u, v and w depend not only on the coordinates x, y 
and z but additionally on the time t, that is they will be func- 
tions of these four quantities x, y, z and t; furthermore, their 
differentials are going to have the following form 

du = Ldx + Idy + Xdz + £dt; 
dv = Mdx + mdy + iidz + DJldt: 
dw = Ndx + ndy + vdz + Vldt; 

Meanwhile we shall always have L + m + v = 0, therefore at 
every arbitrary instant the time t is considered to be constant, 
or dt = 0. Howsoever the functions u, v and w vary with time 
t, it is necessary that at every moment of time the following 
holds: 

du dv dw 
dx dy dz 

Since the condition expresses that any arbitrary portion of the 
fluid is carried in a time dt into a volume equal to itself, the 
same will have to happen, due to the same condition, in the 
next time interval, and therefore in all the following time in- 
tervals. 



the motion of any element is either non-uniform or varying 
in its direction, the change of the motion must be in accor- 
dance with the forces acting on this element. The change of 
the motion becomes known from the known forces, and the 
preceding formulas contain this change; we will now deduce 
new conditions 11 which single out the actual motion among 
all those possible up to this point. 

40. Let us arrange this investigation in two parts as well; 
at first let us consider all motions being performed in the same 
plane. Let AL = x, LI = y be, as before, the defining coor- 
dinates of the position of an arbitrary point I; now, after the 
elapsed time t, the two velocities of the point I parallel to the 
axes AL and AB are u and v: since the variability of time has 
to be taken into account, u and v will be functions of x, y and 
t themselves, in respect of which we put 

du = Ldx + Idy + Zdt andefo = Mdx + mdy + VJldt 

and we have established above that because of the former con- 
dition encoutered above, we have L + m = 0. 

41. After an elapsed small time interval dt the point I is 
carried to p, and it has travelled a distance udt parallel to the 
axis AL, a distance vdt parallel to the other axis AB. Hence, 
to obtain the increments in the velocities u and v of the point I 
which are induced during the time dt, for dx and dy we must 
write the distance udt and vdt, from which will arise these 
true increments of the velocities 

du = Ludt + Ivdt + Zdt and dv = Mudt + mvdt + Tldt. 

Therefore the accelerating forces, which produce these accel- 
erations are 

Accel, force w.r.t. AL = 2(Lu + Iv + £) 
Accel, force w.r.t. AB = 2(Mu + mv + M) 

to which therefore the forces acting upon the particle of water 
ought to be equal. 12 

42. Among the forces which in fact act upon the particles 
of water, the first to be considered is gravity; its effect, how- 
ever, if the plane of motion is horizontal, amounts to nothing. 
Yet if the plane is inclined, the axis AL following the inclina- 
tion, the other being horizontal, gravity generates a constant 
accelerating force parallel to the axis AL, let it be a. Next we 
must not neglect friction, which often hinders the motion of 
water, and not a little. Although its laws have not yet been ex- 
plored sufficiently, nevertheless, following the law of friction 
for solid bodies, probably we shall not wander too far astray 
if we set the friction everywhere proportional to the pressure 
with which the particles of water press upon one another. 13 



II. SECOND PART 

39. Having exposed what pertains to all possible motions, 
let us now investigate the nature of the motion which can re- 
ally occur in the fluid. Here, besides the continuity of the fluid 
and the constancy of its density, we will also have to con- 
sider the forces which act on every element of the fluid. When 



11 Here Euler probably has in mind the condition of potentiality, which he 
will obtain in §§. 47 and 54 for the two-dimensional case and in §. 60 for 
the three-dimensional case. 

12 The unusual factors of 2 in the previous equations have to do with a choice 
of units which soon became obsolete; cf. Truesdell, 1954; Mikahailov, 
1999. 

13 It is actually not clear why Euler takes the friction force proportional to the 
pressure. 
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43. First, must be brought into the calculation the pres- 
sure with which the particles of water everywhere mutually act 
upon each other, by means of which every particle is pressed 
together on all sides by its neighbors; and in so far as this pres- 
sure is not everywhere equal, to that extent motion is commu- 
nicated to that particle. 14 The water simply will be put ev- 
erywhere into a state of compression similar to that which 
still water experiences when stagnating at a certain depth. 
This depth is most conveniently employed for representing the 
pressure at an arbitrary point I of the fluid. Therefore let that 
height, or depth, expressing the state of compression at I, be p, 
a certain function of the coordinates x and y, and should the 
pressure at I vary also with the time, the time will also enter 
into the function p. 

44. Thus let us set dp = Rdx + rdy + Dldt, and let us 
consider a rectangular element of water, Imno, whose sides 
are Im = no = dx and In = mo = dy, whose area is dxdy 
(Fig. 3). The pressure at I is p, the pressure at m is p + Rdx, at 
n it is p+ rdy and at o it is p + Rdx + rdy. Thus the side Im is 
pressed by a force = dx(p + ^Rdx), while the opposite side 
no will be pressed by a force = dx(p+^Rdx+rdy); therefore 
by these two forces the element Imno will be impelled in the 
direction In by a force = —rdxdy. Moreover, in a similar 
manner from the forces dy(p + ^rdy) and dy(p + Rdx + 
\rdy), which act on the sides In and mo will result a force 
= —Rdxdy impelling the element in the direction Im. 




equations: 15 

a - R = 2Lu + 2lv + 2£ or R = a - 2Lu - 2lv - 2£ 
-r= 2Mu + 2mv + 2WI and r = -2Mu - 2mv - 29JI 

from which we gather that 

dp = adx - 2(Lu + lv + Z)dx - 2(Mu + mv + M)dy + VHdt, 

a differential which must be complete or integrable. 

46. Because the term a dx is integrable by itself and noth- 
ing is determined for 9t, from the nature of complete differ- 
entials it is necessary that the following holds in the notation 
already employed: 

d.Lu + Iv + £ _ d.Mu + mv + Wl 
dy dx 

Since ^ =L, % = I; % =M, and % = m it follows that 



, udL , vdl dZ udM vdm dffl 

L/+— — +lm+— +— = MLH — \-mM+ — — = — 

dy dy dy dx dx dx 

which is reduced to this form: 

(L + m)(/-M)+ 

ZdL_rfM\ (^_<bn\ ^_^_ () 
\dy dx J \dy dx J dy dx 

47. In fact, since we knew Ldx + Idy + Stdt and Mdx + 
mdy + Wldt to be complete differentials, 

d~L _ dl dm _ dM^ dZ _ dl Wl _ dM 

dy dx' dx dy ' dy dt dx dt 

after the substitution of which values we have the following 
equation 

(L + m)(/-M)+ 

'dl-dM\ fdl-dM\ dl-dM 
u | : \+v \ ; I H — = 0. 



dx 



dy 



dt 



Plainly, this is satisfied if I = M: so that ^ = ^. Since this 
condition requires that ^ = ^|, 16 it appears finally that the 
differential formula udx + vdy must be complete; in this lies 
the criterion of actual motion. 

48. This criterion is independent from the preceding one, 
which was provided by the continuity of the fluid and its uni- 
form constant density. Therefore even if the fluid in motion 



45. Thus will originate an accelerating force parallel to 
Im = R and an accelerating force parallel to In = —r, 
of which the one directed along the force of gravity a gives 
a — R. Having ignored friction so far, we obtain the following 



14 Here Euler makes full use of the concept of internal pressure, cf. Darrigol 
and Frisch, 2008. 



15 Here the so-called Euler equations of incompressible fluid dynamics appear 
for the first time, but the notation and the units are not very modern, in 
contrast to the memoir he will write three years later (Euler, 1755). 

16 Here there are two problems. The minor problem is a typographical error 
in the published version instead of ^), which is not present in a 
1752 copy of the manuscript (not in Euler's hand), henceforth referred to as 
Euler, 1752. A more serious problem is that Euler here repeats the mistake 
of d'Alembert, 1752 who confused a sufficient condition - the vanishing 
of the vorticity - with a necessary one. 
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changes its density, as happens in the motion of elastic flu- 
ids such as air, this property will hold nonetheless, namely 
udx + vdy has to be a complete differential. In other words, 
the velocities u and v must always be functions of the coordi- 
nates x and y, together with time t, in such a way that when 
the time is taken constant the formula udx + vdy admits an 
integration. 

49. We shall now determine the pressure p itself, which is 
absolutely necessary for perfectly determining the motion of 
the fluid. Since we have found that M = I we have 

dp = adx — 2u(hdx + Idy) — 2v(ldx + mdy) — 2£dx 

-2Wldy + mt. 



Moreover Ldx + Idy 
hence we have 



du — £dt; ldx + mdy = dv — SDtdi; 



dp = adx— 

2udu - 2vdv + 2Zudt + 2Wlvdt - 2Zdx - 2Tldy + mdt 

Therefore, if we wish to ascertain for the present time the pres- 
sure at each point of the fluid, with no account of its variation 
in time, we shall have to consider this equation 

dp = adx — 2udu — 2vdv — 2£,dx — 2Wldy, 



and in our notation £ = ^ and 971 



§ , 17 Hence 



, , du n dv 

dp = adx — 2uau — 2vdv — 2— ax — 2— dy, 
1 dt dt y ' 

in the integration of which the time is to be taken constant. 

50. This equation is integrable by hypothesis, and is in- 
deed understood as such, if we consider the criterion of the 
motion which, as we have seen, consists in that udx + vdy 
be a complete differential when the time t is taken constant. 
Let therefore S be its integral, which consequently will be a 
function of x, y and t themselves. For dt = we obtain 
dS = udx + vdy, while assuming the time t variable as well, 
let us write 



dS = udx + vdy + Vdt, 



on which account we obtain % 



factU = f. 



(ii. 



^ and # 

dx dt 



Then, in 

ay 



51. After inserting these values we will obtain 

du dv , dU , dU , 
-jr-dx + — .dy = — -dx + —.dy 

dt dt dx dy 

and this differential formula is manifestly integrated at con- 
stant time t to give U. For this to become clearer, let us set 

K and ^ = k, so that 



dU 



Kdx 

dXJ 



kdy; thus ^ = ^ — t 

a ' dx dy 

^■dx+^ = Kdx+kdy = dU. Since its integral is U = w 
we shall have 



dS 



from where it appears by integration: 

2dS 

p = Const. + ax — uu — vv — 

dt 

with a given function S of the coordinates x, y and t them- 
selves, whose differential, for dt = is udx + vdy. 

52. In order to understand better the nature of these for- 
mulas, let us consider the true velocity of the point I, which 
is V = \f(uu + vv). And the pressure will be: p = 
Const. + ax - VV - ^ : in which the last term dS de- 
notes the differential of S = J (udx + vdy) itself, where the 
time t is allowed to vary. 

53. If we now wish to also take into account friction, let us 
set it proportional to the pressure p. While the point I travels 
the element ds, the retarding force arising from the friction is 
= j ; so that, setting ^ = U, our differential equation will be 
for constant t 

dp = adx-jds- VdV - 2dU, 

from where we obtain by integration, taking e for the number 
whose hyperbolic 18 logarithm is = 1, 



p = e~f J e l (adx - 2VdV - 2dV) 



or 



p = ax — VV — 2U — - eT J el (ax - VV - 2V)ds. 

54. The criterion of the motion which drives the fluid in 
reality consists in this that, fixing the time t, the differential 
udx + vdy has to be complete: also continuity and constant 

= 0, hence it follows 



uniform density demand that 



dv 

dx 1 dy 

too that this differential udy — vdx will have to be complete. 19 
From where both velocities u and v jointly must be functions 
of the coordinates x and y with the time t in such a way that 
both differential formulas udx + vdy and udy — vdx 20 be 
complete differentials. 

55. Let us set up the same investigation in general, giving 
the point A three velocities directed parallel to the axes AL, 
AB, AC. Let u, v, w denote these functions, which depend on 
coordinates x, y, z, besides t. After a differentiation we obtain 

du = Ijdx + Idy + Xdz + Zdt 
dv = Mdx + mdy + adz + DJtdt 
dw — Ndx + ndy + vdz + Vldt. 

Although here the time t is also taken as variable, nonetheless 
for the motion to be possible, by the preceding condition 21 we 
have L + m + v = 0, or, which reexpresses the same 



du dv dw 
dx dy dz 



= 0, 



dp = adx — 2udu — 2vdv — 2dXJ 



The printed version has L = 2| instead of £ = 2|. Euler, 1752 is 
correct. 



8 Natural 

9 The published version has udx+vdy, a mistake not present in Euler, 1752. 
:0 Previous mistake repeated in the published version. 

: 1 From Part I. 
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a condition on which the present examination does not de- 
pend. 

56. After the passage of time interval dt the point A is 
carried to ir, and it travels a distance udt parallel to the axis 
AL, a distance vdt parallel to the axis AB and a distance wdt 
parallel to the axis AC. Thus the three velocities of the point 
which has moved from A to ir will be: 

parallel to AL = u +Lu dt + l v dt + Xwdt + £ dt; 
parallel to AB — v + Mu dt + mu dt + pw dt + OJldt; 
parallel to AC = w + Nu dt + nv dt + v w dt + 9t dt, 

and the accelerations parallel to the same directions will be 

par. AL = 2(Lw + / v + Aw + £); 
par. AB = 2 (Mi + mv + pw + 9Jt); 
par. AC = 2(Nm + nv +vw + Vl). 

57. If we take the axis AC to be vertical, in such a way that 
the remaining two AL and AB are horizontal, the accelerating 
force due to gravity arises parallel to the axis AC with the 
value —1. Then indeed, denoting the pressure at A by p, its 
differential, at constant time is 

dp = R dx + rdy + pdz, 

from which we obtain the three accelerating forces 

par. AL = R; par. AB = — r; par. AC = — p 

which are in fact easily collected in the same manner as was 
done in §§. 44 and 45, so that it is not necessary to repeat 
the same computation. Hence we obtain the following equa- 
tions 22 

R= -2(Lu + lv + Au> + £) 
r = -2(Mu + mv + pw + 971) 
p = -1 - 2(Nw, + nv + vw + m) 

58. Since the differential formula dp = Rdx + rdy + pdz 
has to be a complete differential, we have 

dR dr dR dp dr dp 

dy dx ' dz dx 7 dz dy 

After a differentiation and a division by —2 the following three 
equations are obtained 23 

*nr + t + tgt + f + u + lm + Xn 

^ + ^ + ^ + lf+ML + mM + uN 

u , ^ + v -i + + £ + la + In + Xu = 
i ^ + i ^ + : lr + ; lF + NL + nM + iyN 



59. Moreover, because of the nature of the complete dif- 
ferentials, we have 



III 



udM _|_ vdm , wdfi . d.VJX 

~T ri-y T" ~T 



dz 
udN 



I vdn I 
dy 1" dy + 



dz 
u>dv 



dz 
dm 



H- MA + mfi + jiv 
+ N/ + rim + vn 



<JL 


dl 


dm, 


dM 


dX 


dl 


flu 


dx ' 


dx 


rln ' 
uy 


nil 
ay 


~ dz ' 


dp 


dM 


d£ 


dl 


dM 


dM 


dx 


dz 


dy 


dt ' 


dx 


dt 


dL 


d\ 


dl 


dX 


dn 


dN 


dz 


dx ' 


dz 


dy ' 


dx 


dy ' 


dv 


dN 


d£ 


dX 


dm. 


dN 


dx 


dz ' 


dz 


dt ' 


dx 


dt ' 


dM 


dp 


dN 


dn 


dm 


dp 


dz 


dx ' 


dy 


dx ' 


dz 


dy ' 


dv 


dn 


dTl 


dp 


dm 


dn 


dy 


dz ' 


dz 


dt ' 


dy 


dt ' 



after substituting of which values those three equations will 
be transformed into these 24 



w 



w 



w 



dl-dM 
Jt 

'dl-dM 
dz 
dX-dN 
Jt 

dX-dN 
dz 
dp — dn 
Jt 

dp — dn 
dz 



dl - dM 



dl - dM 



dx J \ dy 
+ M)(L + to) + An - /iN = 0, 
dX-dN\ fdX-dN^ 



+ u 



dx J \ dy 
+ (A-N)(L + i/)+//i-nM 

(dp — dn\ ( dp — dn 
dx J V dy 

+ (p - n)(m + v) + MX - Nl 



0, 



0. 



60. Now it is manifest that these three equations are satis- 
fied by the following three values 



/ = M; A = N: 



p 



in which is contained the criterion furnished by the consider- 
ation of the forces. Here therefore follows that in the notation 
chosen we have 25 



du dv 
dy dx ' 



du dw 
dz dx ' 



dv 
dz 



dw 
dy 



these conditions moreover are the same as those which are 
required in order that the formula udx + vdy + wdz be a com- 
plete differential. From which this criterion consists in that 
the three velocities u, v and w have to be functions of x, y and 
z together with t in such a manner that for fixed constant time 
the formula udx + vdy + wdz admits an integration. 



These are the three dimensional Euler equations. 

The printed version contains mistakes not present in Euler, 1752: in the 
formula labelled II, instead of L there is £; in the formula labelled III there 
is a v instead of u. 



These are the equations for the vorticity. 

Here Euler repeats the mistake of assuming that the only solution is zero- 
vorticity flow; in Euler, 1755 this will be corrected. 
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61. Taking the time t constant or dt = 0, we have 

du = Ldx + Mdy + Ndz 
dv = Mdx + mdy + ndz 
dw = Ndx + ndy + vdz 

moreover, for R, r and p the values are 

R = -2(Lu + Mv + Nw + £) 
r = -2(Mu + mv + nw + 9Jt) 
p = -1 - 2(Nw, + nv + vw + 91) 

Regarding the pressure p, we obtain the following equation 

dp = —dz 

- 2u(Ldx + Mdy + Ndz) = -dz - 2udu - 2vdv - 2wdw 

- 2v(Mdx + mdy + ndz) - 2Zdx - 2WI - 2 < Jldz 

- 2w(Ndx + ndy + vdz) 

- 2Zdx - 2Tldy - 29Tdz 

62. Since in truth £=^;9Jt=§;9?=^f, we obtain 
by integration 



p = C — z — uu — vv — WW 



/( du , dv , dw , 
{-dt dx + H dy + H dz 



By the previously ascertained condition udx + vdy + wdz 
is integrable. Let us denote its integral by S, which can also 
involve the time t; taking also the time t variable, we have 

dS = udx + vdy + wdz + XJdt, 

and we have f = § = f ; % = £. From where, 
with time generally taken constant, it can be assumed in the 
above integral that 

dU , dU , dU , 

— — dx + — — dy H — — dz — dU, 

dx dy dz 

and we obtain 26 

p = C — z — uu — vv — ww — 2U, or 

p = L — z — uu — vv — ww — 2— 

dt 

63. Thus, uu+vv+ww is manifestly expressing the square 
of the true velocity of the point A, so that, if the true velocity 
of this point is denoted V, the following equation is obtained 
for the pressure 27 

p = C-*-W-^. 

dt 



To use this, firstly one must seek the integral S of the formula 
udx + vdy + wdz which should be complete. This is differen- 
tiated again, taking only the time t as variable. After division 
by dt, one obtains the value of the formula which enters 
into the expression for the state of the pressure p. 

64. But before we may add here the previous criterion, re- 
garding possible motion, the three velocities u, v and w must 
be such functions of the three coordinates x, y and z, and of 
time t that, firstly, udx + vdy + wdz be a complete differen- 
tial and, secondly, that the condition ^ + ^ + ^7=0 holds. 
The whole motion of fluids endowed with invariable density 
is subjected to these two conditions. 

Furthermore, if we take also the time t to be variable, and 
the differential formula udx + vdy + wdz + XJdt is a complete 
differential, the state of the pressure at any point A, expressed 
as an altitude p, will be given by 

p = C — z — uu — vv — ww — 2U, 

if only the fluid enjoys the natural gravity and the plane BAL 
is horizontal. 

65. Suppose we had attributed another direction to the 
gravity or even adopted arbitrary variable forces acting on the 
particles of the fluid. Differences would arise in the values of 
the pressure, but the law which the three velocities of the fluid 
have to obey would not suffer any changes. Thus, whatever 
the acting forces, the three velocities u, v and w have to satisfy 
the conditions that the differential formula udx+vdy+wdz be 
complete and that ^ + ^ + ^ = should hold. Therefore, 
the three velocities u, v and w can be fixed in infinitely many 
ways while satisfing the two conditions; and then it is possible 
to prescribe the pressure at every point of the fluid. 28 

66. However, much more difficult would be the following 
question: given the acting forces and the pressure at all places, 
to determine the motion of the fluid at all points. Indeed, we 
would then have some equations 29 of the form p = C — z — 
uu — vv — ww — 2U, from which the relation of the functions 
u, v and w would have to be defined in such a way that not 
only the equations themselves would be satisfied , but also the 
previously contributed rules 30 would have to be obeyed; this 
work would certainly require the greatest force of calculation. 
It is fitting therefore to inquire in general into the nature of 
functions proper to satisfy both criteria. 

67. Most conveniently therefore let us begin with the char- 
acterization of the integral quantity S, whose differential is 
udx + vdy + wdz, when time is held constant. Let thus S be 
a function of x, y and z, the time t being contained in con- 
stant quantities. When S is differentiated, the coefficients of 
the differentials dx, dy and dz are the velocities u, v and w 
which at the present time suit the point of fluid A, whose co- 
ordinates are x, y and z. The question thus arises here to find 



The published version has a ds in the denominator, instead of the correct 
dt, found in Euler, 1752. 

This is basically the Bernoulli pressure law for potential flow. 



Many statements in this paragraph are rendered invalid by the generally 
incorrect assumption of potential flow. 

The plural is here used probably because this relation has to be satisfied at 
all points. 

Incompressibility and potentiality. 
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the functions S of x, y and z such that + ^ + ^ = 0; 
now, since we have m = w = ^ and to = ^| it follows 
that ^| + + t§ = 0. 31 

da;' 1 dy z dz z 

68. Since it is not plain how this can be handled in general, 
I shall consider certain rather general cases. Let 



S = (Az + By + Cz) n . 



We have 



dx 
ddS 

dx 2 



= nA{Ax + By + Cz)™- 1 and 
= n(n - 1)AA(A# + By + Cz) 



n-2 



and the expressions for ^| and will be similar. Thus we 
have to satisfy 



n(n - l){Ax + By + Cz)"~ 2 (AA + BB + CC) = 

which is plainly satisfied when either n = or n = 1. Thus 
we have the solutions S = Const, and S = Ace + By + Cz, 
where the constants A, B and C are arbitrary. 

69. But if n is neither 0, nor 1, we necessarily have: AA + 
BB + CC = 0: and then S is given by 

S = (Ax + By + Cz)" 

for any value of the exponent n; even the time t itself will 
possibly enter in n. Furthermore we can add up arbitrarily 
many such S and obtain yet another solution. 32 The function 

S = a + /3x + -fy + Sz + e(Ax + By + Cz) n + 
((A'x + B'y + G'zf + V (A"x + B"y + C" z) n " + 
e(A"'x + B"'y + C"'z) n "' etc. 

will satisfy the condition only if we have: 

AA + BB + CC = 0; A' A' + B'B' + C'C = 0; 
A" A" + B"B" + C"C" = etc. 



70. Here suitable values are given for S in which the co- 
ordinates x, y, z have either one, or two, or three, or four 



dimensions 33 
I. S = A 

II. S = Ax + By + Cz 

III. S = Axx + Byy + Czz + 2B>xy + 2Exz + 2Fyz 
with A + B + C = 

IV. S = A.t 3 + By 3 + Cz 3 + 3D:ra;y + 3F:ra;z + Hyyz+ 
6Kxyz + SExyy + ?>Gxzz + Slyzz 

with A + E + G = 0; B + D + I = 0; 
C + F + H = 

V. 

+ Ax 4 + 6Dxxyy + 4Gx 3 y + 4Hxy 3 + lZNxxyz 
S = + By 4 + 6Ea;a;zz + 4Lr 3 z + 4Kxz 3 + YlOxyyz 
+ Cz 4 + 6Fyyzz + 4Ly 3 z + 4Myz 3 + 12Pxyzz 
with A + D + E = G + H + P = 
B+D+F=0 I+K+O=0 
C + E + F = L + M + N = 

71. Hence it is clear how these formulas are to be got- 
ten for any order. First, simply give to the various terms the 
numerical coefficients which belong to them from the law of 
permutation, or, equivalently, which arise when the trinomial 
x + y + z is raised to that same power. Let indefinite letters 
A, B, C, etc., be adjoined to the numerical coefficients. Then, 
ignoring the coefficients, observe whenever there occur three 
terms of the type LZ.t 2 + MZy 2 + NZz 2 having a common 
factor Z formed from the variables. Whenever this occurs, set 
the sum of the literal coefficients L + M + N equal to zero. 
For example, for the fifth power we have 



^IL Lll^ 111 111 \>\J VV V^l VVV^ lltlV^ 

5Da; 4 y + 5Dx i z + lOGx 3 yy + <&x z zz+ 
vz + 30Na;yyzz+ 

)x 4 y + 5<Bx 4 z + 10Bx 3 yy + Sjx 3 zz+ 
20Lx 3 yz + 30Oxyyzz 

+ Gx 5 + 5Fx 4 y + 5$x 4 z + I01x 3 yy + 3x 3 zz+ 
20Mx 3 yz + SOPxyyzz 



S =Aa; 5 + 
20Ka; 3 y 
Bx 5 



and the following determinations of the coefficient letters are 
obtained 

A + G + © = 0; D + H + O = 0; D+I + P = 0; 
B + H + jo = 0; E + G + N = 0; © + F + P = 0; 
K + L + M = 0; 

C + I + S = 0; F + + N = O; ^ + ^ + O = 0. 

In the same way for the sixth order such determinations will 
give 15, for the seventh 21, for the eighth 28 and so on. 

72. In the very first formula S = A the coordinates x, y 
and z are clearly not intertwined. Thus the three velocities u, v 



This is what will later be called Laplace's equation. 

In modern terms, Euler is here using the linear character of the Laplace 33 In modern terms we would say "which are polynomials in x, y, z of degrees 
equation. up to four". 
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and w are equal to zero, and hence this describes a quiet state 
of fluid. Also the pressure at an arbitrary point for different 
times will be able to vary in an arbitrary manner. Indeed A 
is an arbitrary function of time and, for a given time t, the 
pressure at the point A is p = C — _ z Through this 
formula is revealed the state of the fluid, when it is subjected 
at an arbitrary instant to arbitrary forces, which nevertheless 
balance each other, so that no motion in the fluid can arise 
from them: where it happens, if the fluid is enclosed in a vase 
from which it can nowhere escape, it is also compressed by 
suitable forces inside. 

73. Moreover, the second formula S = Ax + By + Cz, 
after differentiation, gives these three velocities to the point A: 



u = A; v = B and w — C. 



Thus simultaneously, all points of the fluid are carried by an 
identical motion in the same direction. From which the whole 
fluid moves in the same manner as a solid body, carried only 
by a forward motion. But at different times the velocities as 
well as the direction of this motion are able to be varied in an 
arbitrary way, depending on what the extrinsic driving forces 
require. Therefore, the pressure at the point A at the time t on 
which A, B, C depend, is 34 p = C- z-AA-BB-CC- 

74. The third formula S = Axx + Byy + Czz + 2Dxy + 
2Exz + 2Yyz, where A + B + C = 0, gives the following 
three velocities 35 of the point A: u = 2Ax + 2Dy + 2Ez; 
v = 2By + 2Dx + 2Fz; w = 2Cz + 2Ex + 2Fy, oru> = 
2Ex + 2Fy - 2(A + B)z. Here, at a given instant, different 
points of the fluid are carried by different motions; moreover, 
in the time development an arbitrary motion of a given point 
is permitted, because A, B, D, E, F can be arbitrary functions 
of the time t. Finally, a much greater variety can take place, if 
more elaborate values are given to the function S. 

75. In the second case the motion of the fluid was cor- 
responding to the forward motion of a solid body in which, 
plainly, at any instant the different parts are carried by a mo- 
tion equal and parallel to itself. In other cases the motion of 
the fluid could be suspected to correspond to solid-body mo- 
tion, either rotational or anomalous. It suffices to put forward 
such a hypothesis - beyond the second case - to find that it 
cannot take place. Indeed, in order to happen, not only would 
it be necessary that the pyramid irQpa would be equal, 36 but 



also similar to the pyramid Xpvo, or that the following holds 



7T$ = 


Xp = 


: dx = ^/(QQ + qq + $$) 


7T/5 = 


Xv = 


dy = y/(RR + rr + pp) 


7T(7 = 


Xo = 


dz = ^/(SS + ss + aa) 


<&p = 


pv = 


^(dx 2 + dy 2 ) = 



VT(Q - R) 2 + (q- rf + (* - pf) 
$cr = po = ^/{dx 2 + dz 2 ) = 
v/((Q - S) 2 + (q - s) 2 + (* - a) 2 ) 
pa = vo = ^jidy 2 + dz 2 ) = 
V~((R-S) 2 + (r- S ) 2 + (p-a) 2 ), 



where we applied the values taken from §. 32. 

76. Then the three latter equations, combined with the 
former, are reduced to these: 

QR+qr + ^p = 0; QS + qs + $a = and RS + rs+pa = 0. 

Moreover, if the values assigned in §. 34 are substituted for 
the letters Q, R, S, q, r, s, $, p, a and the vanishing terms for 
the rests are neglected, the three former will give 

1 = 1 + 2Ldt; I + M = 0; 
1 = 1 + 2mdt; A + N = 0; 
l = l + 2vdt; p + n = 0, 

so that we have L = m = and v = 0, M = -I, N = -A 
and n = —p. 

77. Thus, the three velocities of this point A would have to 
be compared to the condition that the following hold 37 

du — Idy + Xdz; 
dv = — Idx + pdz; 
dw = — Xdx — pdy. 

But the second condition demands a motion of the fluid such 
that I = M, A = N and n = p; hence all the coefficients I, X 
and p vanish; also the velocities u, v and w will take the same 
value everywhere in the fluid. Therefore it is plain that the 
motion of the fluid cannot correspond to solid-body motion 
other than pure translational. 

78. To ascertain the effect of the forces which act from the 
outside upon the fluid, it is first necessary to determine those 
forces 38 which are required for effecting the motion which we 
have assumed to exist in the fluid. These are equivalent to the 
forces which in fact work upon the fluid; furthermore we have 
seen above in §. 56 that three accelerating forces are required, 
which are here repeated. If an element of fluid is conceived 



The printed version, but not Euler, 1752, has a missing BB in the formula. 
In both the printed version and in Euler, 1752, the first velocity component 
is mistakenly denoted by a. 
In volume. 



In the printed version, but not in Euler, 1752, there are several sign mis- 
takes. 

Here, internal forces are meant. 
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here, whose volume, or mass is dxdydz, the moving forces 
required for the motion are 



par. AL 
2dxdydz 
par. AB 
2dxdydz 
par. AC 
2dxdydz 



2dxdydz(hu + Iv + Xw + £) = 
du du du du\ 
dx dy dz dt J 

2dxdydz(Mu + mv + uw + Wl) 



dv 



dv 



dv 



w- 



dv 

dx dy dz ' dt 
2dxdydz(Nu + nv + vw + 91) 



dw dw dw dw 
dx dy dz dt 



so that by triple integration the components of the total forces 
which must act on the whole mass of fluid may be obtained. 

79. But since the second condition requires that udx + 
vdy + wdz be a complete differential, whose integral is S, let 
us put as before, with time allowed to vary, dS = udx + vdy - 
wdz + XJdt. Since ^ = ^; ^ = ^; # 

dy dx ' dz dx 7 dt 

moving forces emerge: 39 



4^ those three 

dx 



par. AL = 2dxdydz 



par. AB = 2dxdydz 



par. AL = 2dxdydz 



udu + vdv + wdw + d\J 
dx 

udu + vdv + wdw + dU 
dy 

udu + vdv + wdw + <iU 
dz 



80. Let us set now uu + vv + ww + 2U — T. The function 
T depends on the coordinates x. y, z; take it at a given instant 
of timet: 40 

dT = Kdx + kdy + ndz. 

The three moving forces of the element dxdydz are 41 

par. AL = Kdxdydz 
par. AB = kdxdydz 
par. AC = ndxdydz 

and by triple integration these formulas ought to be extended 
throughout the mass of the fluid; thus forces equivalent to all 42 
and their directions may be obtained. Truly this discussion is 
for a later investigation, which I shall not deepen here. 

81. Furthermore, the quantity T = uu + vv + ww + 
2U, which is analyzed in this calculation, furnishes a simpler 
formula for expressing the pressure through the height p; we 
have indeed p = C — z — T when the particles of the fluid 
are pressed upon solely by the gravity. But if an arbitrary 



39 There is a misprint in the printed version, w instead of - 

40 There is a misprint: u instead of k. 

41 Here is again a misprint: k instead of k. 

42 The pressure forces. 



particle A is acted upon by three accelerating forces which are 
Q, q and <£>, acting parallel to the directions of the axes AF, 
AB and AC, respectively, after a calculation similar to the 
previous one has been carried out, the pressure will be given 
by 

p = C + J (Qdx + qdy + - T. 

Thus it is plain that the differential Q + qdy + &dz must be 
complete, as otherwise a state of equilibrium, or at least a pos- 
sible one, could not exist. That this condition must be imposed 
on the acting forces Q, q and $ was shown very clearly by the 
most famous Mr. Clairaut. 43 

82. Here are, therefore, the principles of the entire doctrine 
of the motion of fluids, which, even if they at first sight may 
seem insufficiently fruitful, nevertheless embrace almost ev- 
erything treated both in hydrostatics and in hydraulics, so that 
these principles must be regarded as having very broad extent. 
For this to appear more clearly, it is worthwhile to show how 
the precepts learned in hydrostatics and hydraulics follow. 

83. Let us therefore consider first a fluid in a state of rest, 
so that we have u = 0, v = and w = 0; in view of T = 2U, 
the pressure in an arbitrary point A of the fluid is 



P = C + 



J (Qdx + qdy + <5>dz) - 2U. 



Here, U is a function of the time t itself which we take as 
constant. Indeed, we investigate the pressure at a given time; 
the quantity U can be included in the constant C, so that we 
obtain 



P = C + 



J (Qdx + qdy + <frdz) 



where Q, q an $ are the forces acting on the particle of water 
A, parallel to the axes AL, AB and AC. 

84. The pressure p can only depend on the position of 
the point A that is on the coordinates x, y and z; it is thus 
necessary that J (Qdx + qdy + Qdz) be a prescribed function 
of them, which therefore admits integration. Thus it is firstly 
clear that in the manner indicated the fluid cannot be sustained 
in equilibrium, unless the forces acting on each element of the 
fluid are such that the differential formula Qdx + qdy + <&dz 
is complete. Thus, if its integral is denoted P, the pressure at 
A will be p = C + P. Therefore, if the only force present is 
gravity, impelling parallel to the direction CA, we shall have 
p = C — z; hence, if the pressure is fixed at one point A, 
the constant C can be obtained. From which the pressure at a 
given time will be defined completely at all points of the fluid. 

85. However, with time passing, the pressure at a given 
place can change; and this plainly occurs, if variability is as- 
sumed for the forces impelling on the water, whose calcula- 
tion cannot be made from those forces which are assumed to 



' Clairaut, 1743. 
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act on each element of the fluid, 44 but in such a way that they 
keep each other in equilibrium and produce no motion. But 
if, moreover, these forces are not subjected to any change, the 
letter C will indeed denote a constant quantity, not depending 
on time t; and at a given location A we will always find the 
same pressure p = C + P. 

86. It is possible to determine the the extremal shape of a 
fluid in a permanent state, when it is subjected to no forces. 45 
Certainly, at the extreme surface of the fluid at which the fluid 
is left to itself and not contained within the walls of the vase 
in which it is enclosed, the pressure must be zero. Thus we 
shall obtain the following equation: P = const; the shape 
of the external surface of the fluid is then expressed through 
a relation between the three coordinates x, y and z. And if 
for the external circumference held P = E, since C = — E, 
in another arbitrary internal location A the pressure would be 
p — P — E. In this manner, if the particles of the fluid are 
driven by gravity only, and because p = C — z, the following 
will hold at for the external surface z = C; from which the 
external free surface is perceived to be horizontal. 

87. Next, everything which has so far been brought out 
concerning the motion of a fluid through tubes is easily de- 
rived from these principles. The tubes are usually regarded as 
very narrow, or else are assumed to be such that through any 
section normal to the tube the fluid flows across with equal 
motion: from there originates the rule, that the speed of the 
fluid at any place in the tube is reciprocally proportional to 
its amplitude. Let therefore A be an arbitrary point of such 
a tube, of which the figure is expressed by two equations re- 
lating the three coordinates x, y and z, so that thereupon for 
any abscissa x the two remaining coordinates y and z can be 
defined. 

88. Let henceforth the cross section of this tube at A be rr; 
in another fixed location of the tube, where the cross secion 
is //, let the velocity at the present time be b; now after the 
time dt has elapsed, let the velocity become b + e?b , so that b 

is a function of time t, and similarly with . Hence the true 

ff b 

velocity of the fluid at A will be at the present time V = 'H^r ■ 
Since now y and z are obtained from the shape of the tube, we 
have dy = rjdx and dz = 9dx; thus the three velocities of the 
point A in the fluid, parallel to directions AL, AB and AC, are 



the integral of which is denoted = S, we have: 



//b 



1 



//b 



V 



rr 7(1 + rp] + 69) ' rr ^f(l + m + 69) ' 

//» 
w = ■= , 

rr ^(l + r]T] + 99) 

f 4 b b 

and hence, uu + vv + ww = VV = — -j — : and rr is function 
of x itself, thus of the dependent variables y and z. 

89. Since udx+vdy+wdz must be a complete differential, 



dS = U* d f +rm+m = m dxVll + m + w) , 

rr ^/{l + nn + 69) 



rr 



Moreover, dxy/(l+r]ri+99) expresses the element of the tube 
itself; if we denote it by ds, we shall obtain dS — f^ ds ; 

J ' rr 

although b is a function of the time, 46 here we fix the time 
and, furthermore, the quantities s and rr do not depend on 

time but only on the shape of the tube; thus we have S — 

y f ffds 

J rr 

90. Turning now to the pressure p which is found at the 
point of the tube A, the quantity U has to be considered; it 
arises from the differentiation of the quantity S, if the time 
only is considered as variable, so that we have U = Thus, 



since the integral formula J U^t does not involve time t, on 

the one hand we shall have ^ = U = ^ / ^r-, and on the 
other hand it will follow from S. 80 that: 



T 



/ 4 bb M6 f ffds 
r 4 + ~dT 



Therefore, after introducing arbitrary actions of forces Q, q 
and $, the pressure at A will be 

p = C+ {Qdx + qdy + ^dz)-—, — / — — 

J r* dt J rr 

This is that same formula which is commonly written for the 
motion of a fluid through tubes; but now much more widely 
valid, since arbitrary forces acting on the fluid are assumed 
here, while commonly this formula is restricted to gravity 
alone. Meanwhile it is in order to remember that the three 
forces Q, q and $ must be such that the differential formula 
Q dx + q dy + <I> dz be complete, that is, admit integration. 
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